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Abstract 


Experimental activities and computational fluid dynamics (CFD) simulation are presented in this paper for investigating the performance of an 
anode-supported solid oxide fuel cell (SOFC). The goal of this work is to assess a commercial CFD code, Star-CD with es-sofc module, to simulate 
the current-voltage (J-V) characteristics with respect to the experimental data. Compiled with the geometry of cell test housing, a 3D numerical 
model and test conditions were established to analyze the anode-supported cell (ASC) performance including current density and temperature 
distributions, fuel concentration, and fuel utilization. After adjusting parameters in the electrochemical model, the simulation results showed good 
agreements with the experimental data. The results also revealed that the power density increased while the fuel utilization decreased as the fuel 
flow rate increased. 

Based on the results, this modeling work will be implemented to analyze the distributions of fuel and oxidant gases for the SOFC stack, to 


minimize the thermal gradients inside the stack, and to optimize the manifold/flow passage in the future. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


A solid oxide fuel cell (SOFC) is a device that converts 
chemical energy into electrical power through electrochemical 
reactions at elevated temperatures. SOFCs offer a potential for 
high efficiency, fuel flexibility and low air pollution, and the pos- 
sibility of coupling the integrated gasification combined cycle 
(IGCC) to develop a SOFC coal-based central power plant with 
cleanness and higher efficiency. 

SOFCs operating at high temperatures have advantages, such 
as preserving the high oxide-ion conductivity, high kinetic activ- 
ity without a catalyst, carbon monoxide being a fuel rather than 
a poison, and the exhausted heat can be recovered to heat incom- 
ing fuel and air. On the other hand, high operating temperatures 
cause some material problems such as stresses at the interfaces 
due to different thermal expansion coefficients of the SOFC 
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components, and induce performance degradation, resulting in 
instability and a cost increase. 

In the development of the SOFC, planar type designs have 
received much attention because of the higher power densities 
compared with tubular-type designs. A planar type anode- 
supported cell (ASC) with a very thin electrolyte drastically 
decreases the ionic transport resistance and enables operation 
at intermediate temperatures to mitigate the extent of the SOFC 
material problems. 

Many experimental and numerical simulation studies were 
developed to analyze the SOFC cell performance [1—9]. In par- 
ticular, commercially computational fluid dynamics simulation 
codes are widely used to analyze complex phenomena. Autissier 
et al. [10] proved subroutines in the commercial software, FLU- 
ENT, to calculate reaction rates, and distributions of current 
density, gas flow, cell temperature, and fuel concentration in 
ASC cells. In their research, electrochemical behavior is exper- 
imentally fitted on small cells and then applied to complex 
geometries. The studies of Yakabe et al. [11] adopted a com- 
mercial computational fluid dynamics tool, Star-CD, to calculate 
the detailed SOFC characteristics such as fuel concentrations, 
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Nomenclature 

Aj pre-exponential factor 

c mean constant-volume specific heat 
Cp constant-pressure specific heat 

d thickness 

D mass diffusivity (m? s7!) 

E open circuit voltage 

Fact activation energy 


Faraday constant (96,487 C mol— ty 
diffusional energy flux 

diffusional flux 

static enthalpy 

heat of formation 

current density (A m~?) 

exchange current density (A m~?) 
electrochemical reaction rate constant 
molecular weight (g) 

pressure (Pa) 

universal gas constant (8.314 J (mol K)~!) 
resistance 

rate of mass production or consumption 
source term 

temperature (K) 

operating voltage (V) 

gas flow rate (ml min!) 

mass fraction 


Sa eS eee a 


Greek letters 


p density (kg cm7?) 
T stress tensor 
Subscripts 

a anode 

c cathode 

eff effective 

H2 hydrogen 

Oo oxygen 

m mth composition of fuel reactant 
Subscripts 

0 reference 


temperature and current density distributions, etc., through a 
simplified electrochemical model. Thermal stress calculations 
were analyzed by another commercial code, ABAQUS. The 
results showed a co-flow pattern is advantageous to mitigate the 
temperature differences, and hence to reduce internal stresses. In 
Yakabe and Sakurai’s model [12], the electromotive force coun- 
terbalance with the oxygen ions was used to calculate the electric 
current flow in SOFC. The effects of variations of channel to rib 
width ratios were investigated. It was found that the current den- 
sity and contact resistance at the interfaces increased with the 
channel to rib width ratio, as a result of a smaller contact area 
between interconnect and electrode. Recknagle et al. [13] used 
Star-CD to investigate the effects of cell flow configurations, 


Electronic 
Load 


Fig. 1. Schematic diagram of the cell test station. 


including co-flow, cross-flow and counter-flow. Under similar 
conditions in fuel utilizations and average cell temperatures, the 
co-flow case would yield smaller temperature gradient. 

The goal of this work is to verify the accuracy of a commer- 
cial CFD code, Star-CD (Version 3.26), with the es-sofc module, 
in a simulation of the J-V characteristics against the experimen- 
tal data obtained from the cell tests. Based on the geometry of 
the cell test housing and test conditions, a 3D numerical model 
of Star-CD was established to analyze the anode-supported cell 
performance including current density distribution, fuel concen- 
tration distribution, and flow utilization. 

Section 2 of this paper illustrates the SOFC cell perfor- 
mance experiment established at the Institute of Nuclear Energy 
Research (INER). Section 3 describes the simulation method for 
a planar type anode-supported cell including the model geom- 
etry, thermo-fluid model and electrochemical model. Section 4 
presents simulation results and a comparison with experiment 
data. Finally, a brief conclusion is given in Section 5. 


2. Experimental description 


The single-cell test station to investigate the electrochemical 
characteristics of a planar type SOFC cell has been set up at 
INER. The schematic diagram of the cell test station is shown 
in Fig. 1. The test station consists of a gas manifold system, an 
electronic load module and a single-cell test housing enclosed 
in a furnace. Hydrogen is humidified and supplied to the anode 
compartment, and oxygen or air is directed to the cathode side as 
the oxidant gas. A circular cell of 5 cm diameter was set between 
fuel and air distributions. 

The air and hydrogen from the upper side and lower side 
separately flowed into the center of the cell electrodes and dis- 
tributed radically across the cell surfaces. The effective electrode 
area of test cell was about 20 cm?. In the test, a heating rate of 
1°Cmin7! is set. The NiO in the anode was gradually reduced 
to nickel in situ during the heating-up process with a fuel gas 
composition of 10% H2/90% N2. As the temperature reached 
to setting point, the H? concentration was gradually increased 
to 100% at steps of 10% for every 30-min period. The cell per- 
formance was then carried out at varied fuel concentration and 
flow rates. 

Fig. 2 illustrates two measured characteristics: J-V and 
I-power curves for the ASC tests operating at 800°C with gas 
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Fig. 2. Effect of fuel gas flow rate on cell performance. 


flow rates of 300 and 500 ml min~!, respectively. It reveals that 
only slight difference of current density for two different gas 
flow rates in the region of higher cell voltage while the differ- 
ences expands the decrease of cell voltage. For gas flow rate 
of 500 ml min~!, the peak power density was 610 mW cm~?, 
which was about 15% higher than that of 300ml min™! gas 
flow rate. These test data of the /-V characteristics were used to 
benchmark the Star-CD simulation. 


3. Simulation method 


The widely used Star-CD code [11-13] and an advanced 
electrochemical module es-sofc, developed by the Pacific North- 
west National Laboratory (PNNL), were selected to compute the 
thermal-fluid behavior and the complex electrochemical reac- 
tion in a planar type anode-supported cell. The es-sofce module 
is employed to analyze the distributions of fuel and oxidant to 
the cell. Additionally, cell voltages and current densities can be 
subsequently calculated. 


3.1. Geometric model 


Fig. 3 shows a schematic of the geometric model for a SOFC 
test cell. The 3D cylindrical model includes the cathode inter- 
connect, platinum mesh, positive-electrolyte-negative (PEN), 
nickel mesh, and anode interconnect. The gas flow configuration 
including the ribs produced by horizontal and vertical channels 
on the anode and cathode interconnects is also shown in Fig. 3. 
The major geometric data for the input to Star-CD are listed in 
Table 1. 


3.2. Thermo-fluid model 


Star-CD code solves the Navier-Stokes equations for each of 
the computational mesh element in the fluid region. To simplify 


O, inlet 


p. 


= 
connect 


H, inlet 


Fig. 3. Simulation modeling of the cylindrical test rig. 


cell performance calculations, the following assumptions are 
made: (1) gas mixtures are modeled as perfect gases; (2) gas 
flow is steady and laminar everywhere in the cell channels; (3) 
anode and cathode are considered to be isotropic porous media. 
Based on those assumptions, the governing equations for the 
conservation of mass, momentum, enthalpy and species in Star- 
CD code are illustrated as follows. 


Mass equation: 
a (1) 


Momentum equation: 


dou; 
ot 


ð 
+ (pu jui — Tij) = —~— + si (2) 
OX j Xi 
In the above equations, source terms Sm and s; represent 
the sums of the body and other external forces, and tj is the 
stress tensor relative to the velocity gradient. The density of 


Table 1 

Input parameter for Star-CD 

PEN thickness 0.6 mm 
Anode thickness 0.575 mm 
Electrolyte thickness 0.005 mm 
Cathode thickness 0.02 mm 
Anode porosity 0.3 
Cathode porosity 0.3 

Rib thickness 1 mm 
Interconnect thickness 24mm 
Operation temperature 1073 K 
Operation pressure l atm 
Relative humidity (both anode and cathode) 100% 
Cell temperature 800°C 
Anode fuel H2 
Cathode fuel O2 

Fuel inlet temperature (both anode and cathode) 790°C 
Pre-exponential factor, A; 128 
Activation energy, Eact 550 
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fuel mixture is given by the ideal gas law: 


P Ym 
=F (= ze) (3) 


m 


where Y, is the mass fraction of a component with molecular 


weight Mm. 

Enthalpy equation: 

doh ð oP oP Ou; 
— + —(phu; + Fh,j) = j F 

y ag = O O 


(4) 


The specific enthalpy h of the fluid was defined as the sum of 
the thermal and chemical components: 


h = cpT ict + X Ym Hm (5) 


where Hm is the heat of formation of constituent m; cp the 
mean constant-pressure specific heat at temperature T with c 
being the reference specific heat at reference temperature To. 
For solids and constant density fluids, such as liquids, Star-CD 
solves the transport equation for the specific internal energy e: 


e= eT — 3T? +S YmHm (6) 


where ¢ is the constant-volume specific heat. 

Species equation: 

dpYm 
ot 


ð 
F —(pu jYm + Fm,j) = Rm (7) 
Ox j 


where Fp is the diffusional flux component, and R, is the rate 
of mass production or consumption during chemical reaction. 


3.3. Electrochemical model 


In Star-CD, electrochemical reactions and net enthalpy 
generation are calculated by the advanced electrochemical 
module, es-sofc. The calculations include combustion of oxy- 
gen with hydrogen, the gas—water-shift reaction, the capability 
for steam reformation of methane and addition of nitrogen 
leaked into the system. In the computation, a specified cur- 
rent density determines the amount of Oz from the electrolyte, 
which is then combusted with H? from the electrolyte current 
(H2 + (1/2)O2 —> H20). 

The two limiting current densities for the chemical reactions 
can be determined as: 


> _ 2FDe(T) Pp, 


; 8 

i RTd, 8 

i= 4FPDeg(T) in ( Psystem ) (9) 
RTd. Psystem = Po, 


where Derf(T) is the effective diffusion coefficient, Py, and Po, 
the partial pressures of hydrogen and oxygen, F the Faraday 
constant, and d, and de are, respectively, the thicknesses of anode 
and cathode. 


Electrical performance is evaluated from the resulting fuel 
and air compositions. The cell voltage is calculated as: 


Vii) = VNernst = VOhmic = Veutler—Volmer F Veathode F Vanode 
(10) 


In Eq. (10), VNernst is the Nernst or open circuit poten- 
tial, Vohmic the Ohmic potential drop, Vgutler-Volmer the 
Butler—Volmer polarization, Veathode the cathode polarization 
and Vanode is the anode polarization. Each item is determined 
as: 


RT Po , cathode 
Vi = In = 11 
Nernst AF ( Pi, anode ) ( ) 
Vonmic = İR; (12) 
RT > _4 i 
VButler—Volmer = aF sinh Zio (13) 
V, a l 1 i (14) 
= — in i 
cathode 4F ic 
RT i 
Vanode = OF In { 1 — i (15) 


In Eq. (13), the exchange current density i9 corresponds to 
the dynamic electron transfer rate at equilibrium and can be 
expressed as [14]: 


— Eact 
RT 


In Eq. (16), the pre-factor, A; and the activation energy Ect 
are properties specific to the electrode—electrolyte interface in 
the electrochemical model. In present study, the two parameters 
were calibrated to fit one of the experimental J-V curves. 

The entire meshed volume of the model is divided into about 
95,000 cells by preprocessing software Gridgen. The solution to 
the governing equations is performed using Star-CD by dividing 
the model domain into a number of cells as control volumes. In 
the finite volume method, the governing equations are numeri- 
cally integrated over each of these computational cells or control 
volumes. The finite volume method exploits a collocated cell- 
centered variable arrangement, where all dependent variables 
and material properties are stored at the cell center. The average 
value of any quantity within a control volume is given by its 
value at the cell center. The model yields for a set of boundary 
conditions in about 2000 iterations. 


ig = Aj; - exp 


(16) 


4. Results and discussion 


The performance curves are shown in Fig. 4. In this work, 
the experimental results of 500ml min™! gas rate were used 
to speculate the values of A; and Eact in the numerical model. 
In Fig. 4(a), a maximum of 9% average current density dif- 
ference is found between simulation results and experimental 
data obtained from cell tests, which indicates a good agreement 
between simulation and experiment results at cell voltage from 
1 to 0.6 V. Based on this model, another 300 ml min7! gas flow 
rate is simulated with the same parameters of A; and E, (see 
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Operating Temperature= 800°C 
gas flow rate =500mlI/min 


Cell Voltage(V) 
Power Density(W/cm’) 


= simulate results 
O experimental results 


0 0.2 0.4 0.6 0.8 1 
Current Density (Alcm’) 


OperatingTemperature= 800°C 
gas flow rate =300mlI/min 


Cell Voltage(V) 
Power Density(W/cm’) 


simulate results 
experimental results 


0 0.2 0.4 0.6 0.8 1 
Current Density (A/cm?) 


Fig. 4. Comparison between experimental and simulate results. (a) V = 
500 ml min! and (b) V = 300 ml min”! 


Table 1) to verify the reliability of those parameters on dif- 
ferent gas flow rate. Fig. 4(b) displays the maximum current 
density difference is about 10% in the comparison of cell per- 
formance curves between experimental and simulation results 
at 300 ml min~! gas flow rate. From Fig. 4, the adjusting param- 


Fig. 5. Velocity vector distribution in anode-side at Vee =0.7 V. (a) V = 500 ml min 
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Fig. 6. Current density distribution in PEN at different cell voltages. (a) V = 
500 ml min“! and (b) V = 300 ml min™!. 


eters in Star-CD electrochemical model simulations are found 
to have a good agreement in comparison with the experimental 
data. The results also show that the cell performance is improved 
with increased gas flow rate. 

In order to investigate the fluid characteristics of fuel, the 
velocity distribution in anode-side at operating voltage 0.7 V 
with V =500mlmin~! and V = 300mlmin™! are, respec- 
tively, presented in Fig. 5(a and b). Fuel flows into the cell 
from the inlet at the center of interconnect. After gas impinging 
on PEN, the velocity vector demonstrates a radial distribu- 
tion because of an open-outlet flow field and group of pin-fins 
design. As shown in Fig. 5, a higher local velocity is found in 
V = 500 ml min™!. For the open-outlet flow field design, uni- 
form velocity vector can be found in the outlet of the flow field 
for similar long flow routes. 

The operating voltage effects on the local current density 
distributions at V = 500ml min~! and V = 300 ml min“! are 
shown in Fig. 6. Results in Fig. 6 show that current density dis- 
tributions decrease annularly from center to outlet as a result 
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T(°C) 


Veet = 0.6V Veel =0.7V 


Fig. 7. Temperature distribution in PEN at different cell voltages. (a) V = 
500 ml min™! and (b) V = 300 ml min™!. 


of decrease of fuel concentration by electrochemical reaction. 
According to the velocity vector distribution, a slight improve- 
ment of current density is found in the channel area than in the 
rib area, since fuel diffuses into electrodes directly in channel 
areas, but fuel has a longer diffusion-path to get into electrode 
to react in rib areas. By comparing the ordinates in Fig. 6(a and 
b), it indicates that the variations in the local current density are 
more noticeable in the cases of lower operating voltage. The 
electrochemical reaction is expected to be stronger at a lower 
operating voltage and thus a higher operating current density. 

Heat generation from electrochemical reaction of SOFC is 
a significant issue in the development of SOFC. Fig. 7 shows 
temperature distributions in PEN area at different cell voltages. 
The predictions in Fig. 7 disclose that the temperature gradi- 
ent increases gradually with cell voltage decrease because of 
higher electrochemical reaction. In this study, the maximum dif- 
ference in temperature distribution is less than 3 °C. Inspection 
of Fig. 7 indicates that weak temperature peaks distribute on the 
area below channel area as a result of higher electrochemical 
reaction rate with higher fuel gas distributions. It is interesting 
to note that the temperature distributions on gas flow rate are 
contrary to the distribution of current density and fuel gas mass 
fraction because a low temperature fuel and high fuel gas veloc- 
ity (formed by a higher gas flow rate) in the flow-field can carry 
more heat out. In addition, from the results we can find that PEN 
temperature depends on the inlet temperature of cathode gas due 
to a larger specific heat. 

In the study of fuel transport in SOFC, exploring the distribu- 
tions of the fuel gases is important to the design of SOFC flow 
field. To this end, the cell voltage effects on the local distributions 
of oxygen and hydrogen mass fractions at the interface between 
interconnect and electrode are presented in Figs. 8 and 9, respec- 
tively. This can be made plausible by noting that, for a lower 
operating voltage, the electrochemical reaction is faster and thus 
consumes more fuel gases both in anode and cathode. Careful 
inspection of Fig. 8(a) indicates that a slightly higher hydrogen 
mass fraction can be found in channel area. The higher hydrogen 


OO; 


Veet = 0.6V Veen = 0.7V 


0.8000C-01 
0.6000E-01 
0.4000E-01 
0.2000E-01 
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Yu, 


Fig. 8. Mass fraction distribution of hydrogen at anode-side porous media at 
different cell voltages. (a) V = 500 ml min~! and (b) V = 300 ml min”! 


mass fraction could provide more chance to let hydrogen react 
in electrode to lead a better cell performance. In a lower gas flow 
rate condition, the obvious reduction on hydrogen mass fraction 
is displayed in Fig. 8(b). The result reflects the lower cell current 
density distribution as shown in Fig. 6. Similar results could be 
found on oxygen mass fraction in Fig. 9, but the advance on 
the mass fraction in channel here is more noticeable than that in 
hydrogen. 

Fuel utilization rate is an important parameter for the devel- 
opment of SOFC. To examine the gas flow rate effect on fuel 
consumption, Fig. 10 shows the comparison of fuel utilization 
rate between V = 500 ml min~! and V = 300 ml min! where 
the fuel utilization rate is defined as: 


fuelout 


Fuel utilization rate = (1 — ) x 100% (17) 


fuelin 


From Fig. 10 we can find that the fuel utilization rate increases 
with the decrease in cell voltage for a higher electrochemical 
reaction rate. It is interesting to note that the fuel utilization rate 
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Fig. 9. Mass fraction distribution of oxygen at cathode-side porous media with 
difference cell voltage. (a) V = 500 ml min™! and (b) V = 300 ml min™!. 
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Fig. 10. Comparison of fuel utilization rate between V = 500 ml min7! and 
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decreases while gas flow rate increases. With the same electro- 
chemical reaction area, higher fuel gas flow rate could provide 
more fuel to react, but also more chemical products are gener- 
ated, which induces resistance for fuel diffusing into electrodes 
to join the chemical reaction. The mechanism also explains the 
cell performance curves in Fig. 2. The slight current density dif- 
ference at high cell voltage with different fuel gas flow rate is 
found since higher fuel utilization in lower fuel gas flow rate. 


5. Conclusions 


The fluid dynamics simulation code, Star-CD with es-sofc, 
was used to verify the performance curve of a solid oxide fuel 
cell acquired by experimental measurement. Good agreements 
have been demonstrated between the simulations and the experi- 
mental data. The detailed characteristics (include current density 
and fuel concentration distribution, and fuel utilization) of SOFC 
were investigated in this study. What follows are the major find- 
ings: 


1. The cell performance was improved with volume flow rate 
increase; a lower volume flow rate means a lower fuel inlet 
velocity and a lower electrochemical reaction rate. 


2. Higher fuel mass fractions allow the hydrogen to react in the 
electrode better for improved cell performance. 

3. Fuel utilization becomes worse as the fuel flow rate increases, 
while the power density is higher. 


This modeling work has the potential to be further extended 
to analyze the distributions of fuel and oxidant for the SOFC 
stack to investigate the mitigation thermal gradients along 
with temperature prediction for the calculation of thermally 
induced stresses, and manifold/flow passage optimization in the 
future. 
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